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Abstract 

We calculate the Bogolubov coefficients for a metric which describes the snapping of a cosmic string. If 
we insist on a matching condition for all times and a particle interpretation, we find no particle creation. 
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1. Introduction. 

The possibility that cosmic strings may be the root of the mechanism explaining galaxy formation is not 
still ruled out [1]. We need further experimental data on the anisotropies in the cosmic microwave radiation, 
lensing of quasar images, or gravitational radiation stemming from the decay of strings to accept or reject 
this alternative to inflationary quantum fluctuations [2]. If cosmic strings exist, they may give rise to vacuum 
fluctuations which in turn, may result in particle production [3-8]. 

In the presence of a time-like Killing vector, one can define in and out states and calculate the Bogolubov 
coefficients to sec whether particle production actually occurs. Another method would be an approximate 
field theory calculation [9] which has been applied to gravitational particle production during string formation 
in [10] with a positive result. 

In this paper we investigate the particle production during the snapping of a cosmic string and contrast 
our result with that of [10]. The background metric is chosen as the Gleiser-Pullin solution [11], which 
describes the snapping of a cosmic string whose ends expand with the velocity of light. 

In Section 2, we set up our notation and discuss the Klein-Gordon equation for a metric. In Section 
3, we solve the Klein-Gordon equation (for a scalar field) and investigate whether particle production occurs 
via the Bogolubov coefficients. 

In this analysis we have used a classical result which allows the definition of in and out states in Kasner 
spaces [12]. The attempts for the computation of the Bogolubov coefficients are given in Sections 3.1 and 
3.2. 

In Sections 3.1, we obtain a solution in terms of the Hankel functions that has a particle interpretation, 
and the behavior as t — > oo can be interpreted as particle production, but this solution is not bounded near 
the v = hypersurface. In Section 3.2, we obtain a continuous solution in terms of hypergeometric functions 
but there is no particle production. 

We recall that the Bogolubov coefficient method has been used in [4] to compute the particle production 
rates in the case of the instantanous disappearence (or appearence) of a cosmic string along the z-axis. In [4], 
the author uses a smoothed form a the metric and finds an exact solution for the wave equation. The 
in and out states correspond respectively to a spacetime with a cosmic string and to an empty spacetime. 
The matching condition for the solutions of the smoothed wave equation gives nonzero transmission and 
reflection coefficients, and this situation is interpreted as particle production. We note however that in [4] 
the number of particles generated and the total energy (obtained by numerical integration) are proportional 
to (a 3 (l — /3 2 ) 2 ) and (a 4 (l — /3 2 ) 2 ) where a is the the string formation rate. Hence they give meaningful 
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physical results only if the increase in the particle production rate is compansated by a decrease in the 
resulting conical defect. In the case of a metric, a — > oo to start with, hence one should not expect a 
finite total particle number and total energy. 

In the Gleiser-Pullin metric that we start with, the usual in and out states correspond to the interior 
of a flat expanding sphere and the exterior of this sphere with a cosmic string along the z— axis. Thus our 
in and out states (see Eq.3.1.4) should be related to Sahni's results [4] in the limit t — > oo. However we 
cannot obtain a continuous solution with apropriate asymptotic behavior. Even in the case we consider a 
discontinuous solution that matches at the top two orders as t — > oo, and obtain "asymptotic Bogolubov 
coefficients" , these lead to a constant particle production rate, hence to an infinite number of particles. 
We belive that these unsatisfactory results are related to using a metric, since we would have similar 
problems in Sahni's approach, if we use a metric, as discussed in the previous paragraph. 



2. The scalar Klein-Gordon equation. 

We start with the metric given by Gleiser and Pullin 

ds 2 = -Adudv + (u- hv) 2 dtp 2 + (u + hvfdz 2 (2.1) 

where h = 1 for v > and h — (3 2 for v < 0, ip is a periodic coordinate, i.e ip £ [0, 2ir] and the points <p> = 
and <p = 2n are identified, z £ (— oo, +oo), u > 0, u < v for u > and u > —f3 2 v for u < 0. Note that h(v) 
has a jump discontinuity, hence h(v)v is continuous, thus the metric in (2.1) is . We can equivalently 
work with a C^ 00 -* form of the metric by letting h(v) to be function that rises fr om (3 2 smoothly on a small 
interval v £ (—a, a). Then the range of v is defined to be h(v)v > u for u > and h(v)v > —u for u < 0. 

The parameter (3 is a number close to but less than 1, and is related to deficit angle of a cosmic string. 
The Ricci tensor is equal to zero, but one component of the Weyl tensor is proportional to a Dirac delta 
function, which indicates the existence of an impulsive wave. 

The Klein-Gordon equation in these coordinates can be obtained as 

1 F 1 F 



{u-hvf w (u+hvr 
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u — hv u + hv 



F u + \ 



+ ■ 



u — hv u + hv 



F v = 



We look for solutions that have e l ( mi P+ kz ) dependence, with m integer to satisfy the periodicity condition. 
Hence we let 

F^e imv e ifcz F, m integer. 
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Then (2.2) reduces to 



m 2 k 2 1 

uv+ (u-hv) 2 + (u + hv) 2 + 2 



1 1 

+ 



u — hv u + hv 



1 

+ 2 



1 1 



hv u + hv 



F u 

(2.3) 



F v =0. 



This equation can easily be solved in the domains v > and v < 0. The problem arises in matching 
these solutions. We would be primarily interested in continuous solutions with oscillatory behavior as u — > oo 
when v — > 0. In the next section it will be seen that this is not possible: continuous solutions have wrong 
asymptotic behaviour and the solutions with correct asymptotic behaviour can be matched only up to two 
leading orders. 

We can then ask whether we can have generalized solutions, i.e., the ones that can have jump disconti- 
nuities or S function discontinuites. Note that a solution with no 8 function discontinuity is bounded. We 
will show below that if F is bounded, then the u dependence of the jump discontinuity is given by Eq.(2.4). 

Assume that F is bounded on (—a, a) and F \ v -_ a = 4>~{u), F \ v — a = <fi + {u). Since h(v)v and its 
derivative are bounded, integrating (2.3) from —a to a (using integration by parts for the last term) and 
letting a — > 0, we obtain 

(<f>t-4>Z) + -(<f> + -r)=0 (2.4) 
u 

Thus (</>+ — (j>~) is proportional to i. 

If F is allowed to have 8 function (and derivatives of 8 functions) discontinuities, then it can be seen 
that any jump discontinuity can be accomodated. But as the physical meaning of such solutions is not well 
defined, we defer the discussion of such solutions. 

In Section 3.1, we obtain solutions with correct asymptotic behavior that can be matched only up to 
O(j). Thus the discontinuity across v — hypersurface obeys (2.4) asymptotically, hence is bounded as 
t — > oo. 



3. Computation of the Bogolubov coefficients. 
3.1. Solution using Hankel functions. 

We can solve the equation (2.3) for v ^ in the following coordinate system. 

t = u + (3 2 v 



(3.1.1) 

s = u— p v 
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This coordinate system results in incoming and outgoing waves for asymptotic times and it is appropriate 
for particle interpretation [12]. Then the equation reduces to 



Id Id d 2 d 2 
llft~i;ih + dt 2 ~~lte 2 



m 2 k 2 
+ 7^ + 



(3 2 s 2 (3 2 t 2 



= 



with solutions 



Z m (as)Zik(at) for v > 
Z™(as)Zik(at) for v < 



(3.1.2) 

(3.1.3a) 
(3.1.36) 



where a is a constant parameter and the Z's are generic functions obeying Bessel's equation. We take 
solutions 

FM=AH£>{a8)H£\at) 
i^(° ut ) =BH£\as)H%\at) + CH < £ ) (as)H[l ) (at) 
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If we want a continuous solution we need to match the respective solutions at t = s. We will try to 
match the solutions asymptotically, we will see that we can do this only in the leading two orders. 



We obtain 



^e'f (-m+ik) 



ir(m + i ) 



2itT(m- \] 



1 T(ik 1 3 



2itT(ik- i) 



1 - 



i m + 1 



^ + o(t- 2 ) 



i r(i 



h i 3 ^ 



1 + 77T 



i r(f + 1) 



2itr(^ 



TT + 0(*-') 
> J 



1 - 



i r(f + 1) 

2itT(f -\) 



ir+0{t- 2 ) 



or 2 ) 



(3.1.5) 



We define 



(3.1.6) 



and 



Then the equation reduces to 



. r(m + f) r(*fc + f) 

w(m, zfc) = ; + — ff 

r(m-|) r(*fc-|) 



A [1 + ^ W (m, tfc) + 0(t- 2 )] =B 1 + ^u;(£, f ) + 0(i" 2 ) 



+ c 



At the first two orders we obtain 



(3.1.7) 



(3.1.8) 



A = B + C 



Auj(m,ik) = (B-C)u(f,f) 
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(3.1.9) 



Using also the condition 

\B\ 2 - \C\ 2 = 1 (3.1.10) 

we solve |i?| 2 and \C\ 2 as 



|o|2_ \V + U>\ |^,2 _ ,„ 1nl 

^ - |o; + a>|2-|a;-a>|2 |G| " |c + ^ - |c - ^ (3 - LU) 

where w = w(m, ifc), w = 

Using the properties of the V function it can be seen that 

|C| 2 =/? 2 (1-^) 2 . (3.1.12) 

Thus the transmission coefficient is constant, hence its integral cannot be finite, and we cannot obtain a 
finite total number of particles. 

As discussed at the end of Section 1, this result is related to using a c(°) metric. 

3.2. Solution using hypergeometric functions. 

In this section we obtain a solution which is continuous across v = hypcrsurface but asymptotically it 
goes as a power of t. Here we use the following unorthodox coordinate system: 

u = pcosh6» (3.2.1) 

P 2 v =psmh6 (3.2.2) 

We note that, since the Fulling et.al. result [12] is shown for systems with Hankel function solutions, a 
positive result in this coordinate system would not mean particle production. Despite this remark, we will 
use our negative result as a plausibility argument, since it is in line with our previous result. 

If we take solutions of the type 

/ = e mi V" 5 (p, 6) (3.2.3) 

our equation reads 

f)2 Q2 a a2 ok 2 9m 2 

(- P *— 2 + 2coth20— -p- - w - w (co^26 - 1) + -^(coth20 + l))g = (3.2.4) 
where k = ik. An immediate solution is 

, . pV _ ^M^^-iV*-*) ,*[* + »!_ i, "S^A _ i, _ 8; (3. 2 . 5) 
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where x = coth20 , 2 Fi(a, b, c, x) is the hypergeometric function. This solution goes as t to a power in the 
asymptotic region; hence lacks the particle interpretation given in [12]. When v = 0, we see that 9 = and 
this reduces to p s . 

For v < we write the solution as, 

This solution also has the power behavior in the asymptotic region. This solution reduces to p s when v = 0. 
Although we achieved the matching without problem with reasonable behaviour at the asymptotic region, 
we see that we can't fix the two independent functions needed. We know that when v < 0, we get two 
independent solutions for the hypergeometric equation 

mi =2 F\(a, b, c; z) (3.2.7) 

u 2 = z 1 - c 2 F 1 [a - c + 1, b - c + 1, 2 - c; z] (3.2.8) 

The second solution is either identically zero or singular at z equals zero depending on the sign of (1 — c). 
We choose the former case. Here c equals one is not allowed. In this case matching can be achieved when 
the second independent function is put to zero. This shows particle production can be consistently put to 
zero. 



4. Conclusion. 

We find that if we insist on a particle interpretation we can not get a consistent matching of our solutions 
in a specific metric. Only if we neglect terms that are 0(t~ 2 ) and further, we can find the matching 
coefficients which give infinite production. Note also that a possible matching is possible in coordinate 
systems where a particle interpretation can not be given. 

We see that we have to be careful in applying approximate method to solutions for metrics with abrupt 
changes. We agree with Bernard-Duncan [13] who state that we have to be very careful in problems using 
(7(0) metrics. We see that any of the methods used in the smoothed case [4], don't seem to work in our case. 
Our result also shows that the results obtained from approximate field theory calculations [9,10] should be 
taken with a grain of salt. 

Ackowlewdgement: We thank Prof.Dr. Yavuz Nutku for suggesting this problem. This work is partially 
supported by TUBITAK, The Scientific and Technical Council of Turkey under TBAG-CG/1. 
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